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The presence of light sterile neutrinos would strongly modify the energy spectrum of the Tritium 
/3-electrons. We perform an analysis of the KATRIN experiment's sensitivity by scanning almost 
all the allowed region of neutrino mass-squared difference and mixing angles of the 3+1 scenario. 
We consider the effect of the unknown absolute mass scale of active neutrinos on the sensitivity of 
KATRIN to the sterile neutrino mass. We show that after 3 years of data-taking, the KATRIN 
experiment can be sensitive to mixing angles as small as sin^ 29^ ~ 10~^. Particularly we show that 
for small mixing angles, aiir^ 26s ^ 0.1, the KATRIN experiment can gives the strongest limit on 
active-sterile mass-squared difference. 



INTRODUCTION 

There is a consensus that nonzero masses of neutrinos 
and the nontrivial mixing between them is the plausible 
framework to explain the outstanding results of plenty of 
neutrino oscillation experiments. The standard approach 
is to have a three active neutrino scenario, with at least 
two nonzero masses and two reasonably large mixing an- 
gles Q. 

An interesting extension to this standard scenario is 
the existence of extra light sterile neutrino states which 
arose for the first time in the light of LSND experiment [2[ 
showing evidence for Vn-^ Ve oscillation (see also the re- 
cent MiniBooNE data [3| that seems to corroborate). An- 
other evidence of the presence of light sterile neutrinos is 
the so-called reactor neutrino anomaly [J| . This anomaly 
is the departure from unity of the ratio of observed rate 
of events to the predicted rate in very short-baseline re- 
actor neutrino experiments. This departure prompted by 
the re-evaluation of v^ reactor flux that revealed 2.5 % in- 
crease in the flux [J] . Also the Gallium anomaly [5| show 
a deficit of v^, produced by intense radioactive sources, 
such that the ratio of observed rate to the predicted rate 
is 0.86 ± 0.05 [5[. All these anomalies can be under- 
stood by adding one (or more) light sterile neutrino states 
to the scenario with three active neutrinos, with active- 
sterile mass-squared difference ArriggL ^ 0.1 eV^. The 
3-1-1 scenario ^6i, lTJ is the simplest scenario to accommo- 
date the presence of light sterile neutrinos. The model 
is composed of 4 flavor states v^,, ot = e,\x,T and s, that 
are the mixture of mass eigenstates Vi with masses m^, 
1 = 1,2, 3, 4. We associate the mass scale that induce the 
very short-baseline oscillations, Ar?T,gy^y , with the mass 
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It should be noticed that 



in order to explain the reactor anomaly, a nonzero com- 
ponent of 1^4 in the state I'e is necessary. 



A new generation of tritium beta decay experiments, 
like the KArlsruhe TRItium Neutrino (KATRIN) ex- 
periment [8[, was proposed to search kinematically for 
the neutrino mass by measuring the energy spectrum 
of the electrons from the beta decay of tritium '^H— >• 
•^He"*" -I- e~ -|- i/e (see also other proposals [9|). A nonzero 
neutrino mass results in displacement of the cndpoint 
energy in electron spectrum which is the focus region 
to probe in KATRIN. In this Letter we analyze the ca- 
pability of KATRIN experiment in the search for end- 
point displacement in the energy spectrum of /3-electrons 
due to the lightest neutrino mass TOi and the irregulari- 
ties in the shape of energy spectrum due to the heavier 
(mostly sterile) neutrino mass m^ and its nonzero mixing 
with electron neutrino. We discuss for the first time the 
role of the lightest neutrino mass r?T,i in determination 
of the sensitivity of KATRIN experiment to the oscilla- 
tion parameter Am^^. We show that with the present 
KATRIN design, this experiment is the only one sensi- 
tive to the part of parameter space corresponding to very 
small active-sterile mixing angle and large active-sterile 
mass-squared difference. 



TRITIUM BETA DECAY AT KATRIN 

The KATRIN experiment [8[ have the following setup. 
Injected molecular tritium gas at Tritium Laboratory 
Karlsruhe provides high luminosity /3-electrons emitting 
isotropically. The electrons will be guided by the gra- 
dient of a magnetic field to the so-called MAC-E-Filter 
(Magnetic Adiabatic Collimation combined with an Elec- 
trostatic Filter) spectrometer. The ratio of the minimum 
magnetic field at the central plane of the spectrometer 
(Ba = 3 X 10~^ T) to the maximum magnetic field near 
the tritium source (-Bmax = 6 T) determines the rela- 



tive sharpness of the energy filtering of MAC-E-Filter. 
Also applying a magnetic field Bs = 3.6 T at the tritium 
source suppresses the entrance of electrons with large ini- 
tial emission angle. In the spectrometer, with the help 
of an electric field parallel to the electron's propagation 
path, it is possible to make an electrostatic barrier qU 
which can be passed just by electrons with energy higher 
than the height of the barrier. Taking all together, the 
transmission function of the KATRIN spectrometer as 
a function of electron kinetic energy Ke and retarding 
potential qU is 
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where AKe/Ke = Ba/ Bra&x is the relative sharpness of 
the filter. However, electrons can undergo inelastic scat- 
tering with tritium molecules in the source which can 
change the spectrum. Taking into account the proba- 
bility of multiple inelastic scattering, the transmission 
function modifies to the following convoluted form 



T'{K,,qU) 



K^-qU 



T{K,-e,qU) 



(1) 



X [PoS{e) + Pif{e) + P2U ® /)(e) + . . .] de, 

where P„ is the probability that the electron scatters n 
times off the tritium molecules before leaving the source 
and /(e) is the energy loss function at each scatter- 
ing [10[ . The symbol ® defines the following convolution: 



(/®/)(e) 



K^-qU 



f{e')f{e-e')de'. 



The rate of the electrons passing the potential barrier 
qU and arriving at the detector is 

S{Q,qU,[U,,],[m,]) = 

J^ /3{K„ g, [Ue^], [m,])T'{K,, qU) dK, , (2) 

where the symbols [rrii,] = {toi, . . . ,to„} and [Uei] = 
{C/ei, • • • , Uen} dcuotes respectively the set of the masses 
of neutrino mass eigenstates Vi and the elements of the 
first row of Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
C^PMNS Hi rnixing matrix. The function /3 gives the 
spectrum of electrons in beta decay 



/3(ife,g,[t/«],K]) = 

N,F{Z, K,)EeP. 






lpifi|C/ej 



£f - m^e(£j - nij) 



where £i = Q — Wi — K^- In the above equation Ef^ and 
Pe are respectively the electron's energy and momentum; 
F{Z,Ke) is the Fermi function which takes into account 
the electrostatic interaction of the emitted electron with 
the daughter nucleus with Z = 2 [12|; Wi and pi are 



respectively the excitation energy and transition proba- 
bility for the excited state i of the daughter nucleus [l3[ ; 
and the Heaviside step function 6 guarantees the con- 
servation of energy. The index i runs over the excited 
states and index j runs over the neutrino mass eigen- 
states. The factor Ng determines the total number of 
emitted electrons which for the KATRIN design param- 
eters is 1.47 X 10-13 g-igV-5 d. 



SENSITIVITY OF KATRIN TO STERILE 
NEUTRINO 

The functional form of the /3-clcctron spectrum 
in Eq. ^ depends on the set of masses [m^l = 



141 , for the case that 



{mi, . . . , TO„}. As discussed in 
the energy resolution of the experiment near the end- 
point and the energy interval that is probed by the ex- 
periment is much larger than the mass splittings, it is 
possible to replace the set of masses [m,y] with an effec- 
tive mass m/3 = \f^~rn^^^J^ . However, in the case of 
sterile neutrino with a mass-squared difference ^ 1 eV , 
this approximation fails and the error of using effective 
mass in the fit of spectrum becomes large. Here we use 
the exact form of Eq. ([3]) with four mass parameters 
{7711,7712,7713,7714} in the case of 3-1-1 scheme. However, 
from these four masses just m\ and 7714 enter the anal- 
ysis and the other two are fixed by TO2 = \pr^^^-Am^^ 
and 7^73 = y^Tif+Aml^; where the values of Am^^ and 
A777§i are fixed by oscillation phenomenology |l|. For 
the elements of the PMNS mixing matrix we use such 
parameterization of the L'4x4 that its 3 x 3 sub- matrix 
for the light active masses reduces to the PDG param- 
eterization [l5[. Thus, the element [/e4 just depends on 
one mixing angle Os^ through |f/e4| = sin 0s . With the 
above mentioned considerations, the rate of the events in 
Eq. ([2]) is a function of parameters {Q,qU,Ue4,mi,m4^). 
The total rate is the sum of signal rate S in Eq. ^ and 
the expected rate of the background events Ni, which 
for the KATRIN is 10 mHz Q. For the Q-value of the 
tritium beta decay we use the central value of a recent 
measurement Q = 18571.8 ± 1.2 eV ^. To illustrate 
the behavior of sterile admixture, we define the following 
ratio 



S'(sin^ 29s,mi,Amli) + Nb 
S'(sin2 261,, = 0, TTii, A7772^ = 0) + Nb 



(4) 



(3) where sin^ 26's = 4|C/e4p(l - |C/e4p). Wc plotted this ra- 
tio in Fig. [T]for different set of the parameters sin^ 29s, 



mi and Amf^. Comparing the black (dotted) curve with 
the red (dashed) curve in Fig. [TJ it is easy to see the 
change in /3-spectrum for different values of the mass 7714 
for a vanishing light mass ttii =0. The height of min- 
imum depends on the values of sin 26s and Amli and 
the position of the minimum depends only on 7774 . How- 
ever, in the comparison between the red (dashed) and 
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FIG. 1. The ratios of the total rate in Eq. Q for various 
mixing and mass parameters. The vertical line shows the 
Q = 18571.8 eV. 



blue (dot-dashed) curves in Fig. [T] we see that by the 
inclusion of a nonzero value for mi the two curves with 
the same Amli cross each other, implying that the lack 
of knowledge about the value of toi can lay a shadow on 



the determination of Ar 



Ml- 



To quantify the sensitivity 



of the KATRIN experiment to the sterile neutrino mass 
we define the following x^ function 



X^{Q,Ue4,rni,m4,Rs,Rb) 



(5) 



E 



(iVexp([gC/]») - iVth(Q, [qU]^, Ue4, mi, mi, Rs,Rb)f 



where a ~ ^/Nc^ is the statistical standard deviation. 
In Eq. ^ the A'th corresponds to the number of de- 
tected /3-electrons when the retarding potential has the 
value [qU]i which is calculated by multiplying the rate in 
Eq. ([2]) by the time spent for the measurement: 

7Vth(. . . , [qU],, . . .) = t[i] ■ {R,S{. . . , [qU]„ ...) + RbNb) , 

where Rs and Rb are respectively the normalization fac- 
tors of signal and background events; and Nb = 10 mHz. 
The iVexp denotes the experimental number of events as- 
suming 7714 and [/e4 equal to zero. The index i in Eq. ([5|) 
runs in 31 steps such that the retarding potential cov- 
ers the range qU e [Q - 20 eV,Q + 5 eV]. For the 
time t[i] spent at each step of the retarding potential 
we use the optimized measurement time proposed by the 
KATRIN collaboration (see Figure 131 in Ref. [1]). We 
minimize the x^ function with respect to the normal- 
ization factors Rs and Rb analytically and with respect 
to Q numerically in its uncertainty range. Fig. [2] shows 
the 90% C.L. sensitivity contours of the KATRIN exper- 
iment in the {sin^ 26 s , Am,fi) plane for the total mea- 
surement time X^i ^W = 3 years. The black (dotted), 
red (dashed) and blue (dot-dashed) curves correspond 
respectively to nii = 0,1,2 cV. The green (solid) curves 
show the 90% C.L. allowed region from the global fit of 
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FIG. 2. The 90% C.L. contours of the KATRIN experi- 
ment in the (sin 26s,Am4i) plane. The black (dotted), red 
(dashed) and blue (dot-dashed) curves correspond respec- 
tively to TTii — 0, 1,2 eV. The green (solid) curves show the 
90% C.L. allowed region and the red cross the best fit point 
for the global fit of the data for the 3-1-1 scheme [31. The 
magenta and purple curves show respectively the BugeyS and 
Bugey4-l-Rovno exclusion curves ^. 



the short-baseline oscillation data [7[ and the red cross 
shows the best-fit value. 

We can conclude from Fig. [5] that after three years of 
data-taking, the KATRIN experiment can exclude the 
main part of the current allowed region of 3-1-1 scenario. 
Beside that, the KATRIN is the most sensitive experi- 
ment in the upper left part of (sin 26s, Anili) plane. As 
can be seen, the experiment is sensitive to mixing angles 
as small as sin^ 26's ~ 10^^. For comparison we have 
shown in Fig. [2] the present exclusion curves of Bugey3 
(magenta solid curve) and Bugey4-|-Rovno (purple solid 
curve) [5|. 

We have also found that, for the first time, varying 
the value of lightest neutrino mass mi can affect the sen- 
sitivity to the large mass 7714. For example, for a fixed 
value of mixing angle sin^ 26s = 0.1, the sensitivity of the 
experiment is Amli — 0.98, 1.1 and 1.5 eV^ for respec- 
tively TOi =0,1 and 2 eV as can be seen from Fig. [2l 
This implies a correlation between the discovery poten- 
tial of Am,\^ and the value of mi. For smaller (larger) 
values of ttii the correlation is weaker (stronger), and 
also the correlation depends on the mixing angle. For 
very small mixing angles the correlation disappear be- 
cause of the weak 777,4 contribution, and for sin 26s = 1 
case, which corresponds to equal admixture of t^i and 1^4 
in Vq, the correlation still exists. Also, it should be no- 



ticed that the large values of mi are in potential conflict 
with the standard ACDM model of cosmology, although 
extensions to non-minimal cosmological models can re- 
duce the conflict |17| . 
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The effect of light sterile neutrinos in beta decay exper- 
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imcnts was formerly discussed in Refs. 
results are in agreement with the estimations of 
In Ref. [19| a general analysis was not performed, but 
we can conclude that we have similar results for small 
values of mi. In Ref. [20| the authors assume null value 
for the lightest mass mi which can be compared with 
the black (dotted) curve in Fig. [2] For mixing an- 
gles sin^ 26's ~ 0.1, our exclusion is ^ 0.6 stronger in 
logio ^'^41 • This stronger exclusion can be the result of 
two issues: i) we use the optimized running time in the 
measurement of spectrum; ii) a different x^ function can 
be used in Ref. [20|. For the very small mixing angles 
sin^ 26s ^ 0.05, our analysis do not show the wiggling 
behavior in Ref. [20| which is not expected in kinemati- 
cal mass measurement experiments. Generally, we have 
an agreement in the limiting case of very small light mass 
mi with the previous results. However, for non small mi 
masses, we found an interplay between the mixing pa- 
rameter C/e4 and the two free mass scales mi and 7714, 
that was not noticed in all previous analyses. 



CONCLUSIONS 

In this work we analyzed the sensitivity of beta de- 
cay experiment KATRIN in determining the mass scale 
associated with the presence of a light sterile neutrino 
state. Motivation comes from the Pe reactor anomaly, 
the Gallium anomaly, the LSND and MiniBooNE exper- 
iments which favor the presence of light sterile neutri- 
nos that mix with electron neutrinos, compatible with 
AmggL ^ 0.1 eV^ and small mixing angles sin^ 26s- 

For the first time we considered the effect of light mass 
scale TTii in the determination of oscillation parameter 
Am^li in KATRIN. We exploited a general treatment of 
nonzero values for the lightest mass scale mi and the 
heavier mass scale 7714. We have shown that varying the 
unknown mass scale ?7ii G [0, 2] eV, induce 0.2 uncer- 
tainty in the sensitivity of KATRIN to log^o Am^j^ . How- 
ever, we have shown that despite this uncertainty, with 3 
years of data-taking, KATRIN can exclude the main part 
of the current allowed region in (sin^ 20^, Am|]^) plane 
indicated by the global fit of short-baseline oscillation 
experiments. Also, we have shown that for very small 
mixing angles sm^26s < 10~^, the KATRIN experiment 
gives the strongest bound on the oscillation parameter 
Aml^. 
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